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We study the charge-density wave phase in TiSe 2 by using first principle calculations. We show 
that, regardless of the local functional used and as long as the cell parameters are in agreement 
with the experiment, density-functional calculations are able to reproduce not only the structural 
instability of TiSe 2 , but also the effective distortion observed in the experiments. We study the elec¬ 
tronic structure evolution of the system under the charge-density wave deformation. In particular, 
we show that the energy bands for the distorted superstructure, unfolded into the original Brillouin 
zone, are in reasonable agreement with angle-resolved photoemission spectroscopy (ARPES) data 
taken at low temperature. On the contrary, the energy bands for the undistorted structure are not 
in good agreement with ARPES at high temperature. Motivated by these results, we investigate 
the effect of the correlation on the electrons of the localized Ti -d orbitals by using the LDA+U 
method. We show that within this approximation the electronic bands for both the undistorted and 
distorted structure are in very good agreement with ARPES. On the other hand, the U eliminates 
the phonon instability of the system. Some possible explanations for this counter intuitive result 
are proposed. Particularly, the possibility of taking into account the dependence of the parameter 
U from the atomic positions is suggested. 


I. INTRODUCTION 


The group IV6 transition metal diselenide lT-TiSe 2 
(space group P3ml) is a layered compound which has 
received considerable attention because of its interesting 
physical properties. In particular, below a critical tem¬ 
perature Xcdw — 200 K it undergoes a commensurate 
charge density wave (CDW) transition, with the forma¬ 
tion of a 2 x 2 x 2 superlattice structure (space group 
P3cl) accompanied by the softening of a zone boundary 
phonon and changes in the transport properties 1 3 . In 
spite of many experimental and theoretical studies, the 
driving force of this structural phase transition remains 
controversial. Several mechanisms have been proposed 
for the origin of the instability in TiSe 2 and they can be 
roughly classified into two main groups depending on the 
driving role played either by the electrons or by the lat¬ 
tice. In fact, a charge density wave occurs always simul¬ 
taneously with a periodic lattice distortions, so with both 
a modification of the electron and phonon spectra, but it 
is unclear if what is observed is primarily an instability of 
the electronic system or of the lattice^. To the first case 
belongs the excitonic insulator modeP®, where the CDW 
is essentially view as a many-body effect originated by the 
poorly screened hole-electron Coulomb interaction giving 
rise to a condensate of excitons and a consequent distor¬ 
tion. In the second fa mily w e find Peierls and Jahn-Teller 
band-type mechanism d 9 * 10 *, where the instability essen¬ 
tially comes from the electron-phonon coupling leading 
the a lattice distortion which lowers the total energy of 
the system. The CDW phase competes with supercon¬ 
ductivity since TiSe 2 is not superconducting at low tem¬ 
perature, but CDW is suppressed and superconductivity 
stabilized either by Cu intercalation 11 or pressure 12 . For 


this reason, a deep and definitive understanding of the 
CDW occurrence would be interesting both for concep¬ 
tual reasons and technological applications. 

In this paper we present the specific role of the 
electron-phonon coupling in the appearance of the CDW 
ordering. As already shown in Ref. [13] thorough density 
functional theory (DFT) calculations it is possible to ob¬ 
serve a structural instability at the L and M points of 
the Brillouin zone (BZ) consistent with a 2 x 2 x 2 (L) and 
2x2x1 (M) real space superstructure. Here we want 
to provide a deeper analysis of this instability in order to 
find the exact distortion predicted and compare it with 
the experimental findings. 

The paper is organized as follows: in Sec. |TT| we sum¬ 
marize the method and the parameters used. Then, in 
Sec.|III[ we present the ab-initio structural analysis of the 
TiSe 2 instability and CDW phase. Afterwards, in Sec.|TV| 
we analyze the electronic structure of the distorted phase. 
Particularly, we unfold the bands of the distorted phase 
into the undistorted BZ in order to compare the theoreti¬ 
cal results with the data from an ARPES experiment. In 
Sec. IVBj we show the results of the LDA+U calculation. 


Finally, conclusions are presented in Sec. VI 


II. COMPUTATIONAL DETAILS 

All calculations were performed within the frame¬ 
work of DFT using the Quantum ESPRESSO pack- 
ag^HI which uses a plane-wave basis set to describe the 
valence-electron wave function and charge density. For 
the exchange-correlation functional we used both the 
Perdew-Zunger local density approximation (LDA) 1 ®and 
the Predew-Burke-Ernzerhof conjugate gradient approx- 
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TABLE I. Experimental and theoretical geometrical parame¬ 
ters of the system in the undistorted phase (cf. Fig0: hexag¬ 
onal lattice constant a, distance c between the layers, distance 
h between the Se and the Ti planes in a layer and horizontal 
projection R of distance between Se and Ti in an octahedron. 
The subscripts ‘Exp’ and c Th’ refer to the experimental and 
theoretical cell, respectively. Notice that for the undistorted 
phase R is fixed by the unit cell geometry (it must be equal to 
ay/ 3/3 in order to obtain a null force along the planar direc¬ 
tion) wh erea s this is not true anymore in the distorted phase 
(cf. Tab.|vH. 



a (A) 

c (A) 

h (A) 

R(k) 

Exp 

3.540 

6.007 

1.532 

2.044 

LDAexp 

3.540 

6.007 

1.499 

2.044 

gga Exp 

3.540 

6.007 

1.534 

2.044 

GGAgw 

3.540 

6.007 

1.532 

2.044 

LDATh 

3.434 

5.792 

1.535 

1.982 

GGA T h 

3.536 

6.719 

1.548 

2.041 

GGAxh W 

3.510 

6.165 

1.553 

2.026 


imation (GGA) 16 . In the second case, as adjacent lay- 
ers in TiSe 2 are coupled by Van der Waals forces, we 
also considered a correction to the functional 1 ^ which is 
aimed at describing more accurately this kind of interac¬ 
tion (GGA VdW ). In these cases the phonons have been 
calculated using density functional perturbation theory 
(DFPT) in the linear responsd^. 

We used a cutoff of 85 Ry and 850 Ry (1 Ry ~ 13.6 eV) 
for the wave functions and the charge density, respec¬ 
tively; the BZ integration has been performed with a 
Monkhorst-Pack gricP^of 24 x 24 x 12 k and a Hermite- 
Gaussian smearing of 0.01 Ry. The self-consistent solu¬ 
tion of the Kohn-Sham equations was obtained when the 
total energy changed by less than 10“ 10 Ry. We studied 
the system with internal theoretical coordinates ( i.e. zero 
theoretical forces) and with both experimental and the¬ 
oretical cell (i.e. zero theoretical pressure). The theo¬ 
retical parameters have been obtained by relaxing the 
structure starting from the experimental parameters 20 
until the forces on the atoms were less than 10 -3 Ry 1 
(ao ~ 0.529177A is the Bohr radius) and the pressure less 
than 0.5 Kbar. The values of the geometrical parameters 
obtained for different local functionals are reported in 
Tab. [Vj| In particular, as it can be seen, the distance be¬ 
tween the layers is underestimated in LDA and overesti¬ 
mated in GGA whereas, as expected, the best agreement 
between theory and experiment is obtained by using the 
GGA vdw functional. 

In order to take into account the strong correlation ef¬ 
fects due to the localized d orbitals of Ti, we also consid¬ 
ered the LDA -\-U method in the simplified form described 
in Ref. 1211221 Since we consider one strong-correlated 
orbital and no spin, this adds a single additional param¬ 
eter U. 

We essentially performed two kind of analysis. On one 
hand we calculated the variation of the phonon frequen¬ 
cies in L and M with the value of U for the experimental 


cell and theoretical internal coordinates obtained by re¬ 
laxing the atomic positions for each value of U. In this 
case the phonon frequencies have been obtained by using 
the finite difference method in Ref. [23l for a 2 x 2 x 2 
cell with al/66~0.015A displacement of the atoms. A 
mesh grid of 24 x 24 x 12 k for the super-cell Brillouin 
zone and a Hermite-Gaussian smearing of 0.125 • 10 -2 Ry 
have been used for the related self-consistent calculations 
of the forces. 

On the other hand we calculated, by linear response, 
a first-principle estimate of U for the LDAexp case in 
the undistorted fase, through the difference between the 
screened and bare second derivative of the total energy 
with respect to the occupation of the Ti -d orbitaP^. For 
an input value U- in , used to define the starting system, 
the linear-response calculation returns a different out¬ 
put value C/out A km but, in order to be consistent 
and replace the LDA interaction term with the corre¬ 
sponding Hubbard correction, the ideal case in which 
U[ n = 7/out = 7/ should be considered. In fact, even if 
it is a common practice to simply compute the ab-initio 
value of U in one step with U- in = 0 , this consideration 
can be relevant, especially if the LDA and LDA +U sys¬ 
tems are qualitative different!^. In our case, in particu¬ 
lar, the effect of U is to open a gap between the bands 
with the result of obtaining, for U ~ 4, a metal-insulator 
transition (see Sec. El For this reason we determined 
U with a self-consistent procedure starting from the un¬ 
perturbed system (U[ n = 0) and using, step by step, the 
obtained U out as U\ n for the subsequent calculation. For 
each step we obtained the result first by performing the 
linear response calculation on a 2 x 2 x 2 cell and then 
extrapolating the outcome to a 6 x 6 x 6 cell (see Ref. 
21 for details). With this procedure we converged in a 
few steps to the value U ~ 3.902 eV. Since we decided to 
work with a fixed configuration, for all the steps we al¬ 
ways kept fixed the internal coordinates equal to the ones 
obtained by relaxing the atomic positions with U = 0; 
moreover, in order to achieve a precision of 10 “ 3 eV for 
the converged value of 77, we set the energy convergence 
threshold for self-consistency equal to 10 “ 14 Ry. 

III. STRUCTURAL ANALYSIS OF THE CDW 

In the BZ of TiSe 2 the three L points - Li, L 2 , £3 
- and the three M points - Mi, M 2 , M 3 - are equiv¬ 
alent thanks to the three-fold rotation symmetry of the 
system (see Fig. [2|. The vectors q L . and q M . from T 
to Li and M^, respectively, have reduced components 
(cfr. Fig. [ 2 ]): 
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FIG. 1. (color online) Octahedral structure of TiSe 2 in a layer. 
The Se atoms in the upper and lower plane are on circles of 
radius R at distance h from the Ti atoms plane. The value 
of the hexagonal lattice parameter a, the distance c between 
two layers (not shown) and h completely define the system 
(the value of R is fixed by the geometry in the undistorted 
phase). The rotational displacement of the Se atoms in a 
triple-g mode and the attraction (repulsion) exerted over the 
Ti atoms by two close Se’s for v > 0 (y < 0) is also shown 
(see the main text for the definition of v). 


In Li and Mi the small group is C^h and the decom¬ 
position in irreducible representations is: 

2 A u © 2 A g © 4 B u © B g (2) 

By using the Density Functional Perturbation Theory 
(DFPT) we computed the phonon frequencies at M and 
L. We found that in all the cases the lowest phonon mode 
has symmetry A u with frequency always imaginary ex¬ 
cept in the LDA T h case. The values for the phonon fre¬ 
quencies in the GGA^h W case are shown in table Tabjn] 
and the frequencies of the lowest mode for all the cases 
are reported in table Tab|ml When a phonon frequency 
l jj is imaginary we conventionally indicate it with the neg¬ 
ative value — |oj | . 


A. Ab-initio analysis of the structural instability 

The imaginary phonon frequencies in M and L corre¬ 
spond to a structural instability consistent with a 2 x 2 x 1 
(M) or a 2 x 2 x 2 ( L ) real-space superstructure. In or¬ 
der to study the 2x2x2 distortions which lower the 
energy we consider the corresponding supercell of the 
undistorted crystal (which has 24 atoms) and the 72 di¬ 
mensional space V whose general element d = di a is the 
displacement of the i-th atom of the supercell along the 
cartesian coordinate a. 



FIG. 2. Brillouin zone of TiSe 2 . The right hand figure is the 
BZ as seen from above. The bi are the reciprocal lattice basis 
vectors. 


TABLE II. Phonon frequencies with relative mode symmetries 
at L and M in the the GGA^h W case. 


ujl (meV) 

c om (meV) 

A u 

-10.00 

-9.33 

B u 

12.18 

12.69 

A u 

14.27 

14.62 

B g 

17.37 

17.25 

B u 

17.95 

17.57 

a 9 

20.33 

20.44 

Ag 

23.67 

23.48 

B u 

24.78 

24.93 

B u 

37.00 

37.84 


TABLE III. Phonon frequencies of the lowest mode ( A u ) at 
L and M for the cases analyzed. 


ujl (meV) 

ujm (meV) 

LDAexp 

-10.38 

-9.17 

gga Bxp 

-9.83 

-8.28 

GGAgw 

-9.61 

-8.03 

LDAih 

+4.32 

+7.13 

GGAxh 

-13.14 

-13.18 

GGAxh W 

-10.01 

-9.34 
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For the 2 x 2 x 2 superstructure the eight points T, A, Li 
and Mi of BZ all refold to the T point. Thus, the space V 
is equal to the orthogonal sum of the corresponding nine 
dimensional subspaces Vr, Va, Vl. , VmV 


V = Vr 0 Va 0 Vl 0 Vm (3) 


3 3 

U = 0 V Li V M = 0 V Mi (4) 

*=1 i =1 

whose vectors describe distortions with a definite mod¬ 
ulation character with respect to the original lxlxl 
unit cell of the undistorted phase. In particular Vi,. and 
Vmi are made of plane-wave lattice distortions with wave 
vector q L . and q M , respectively, that is distortions hav¬ 
ing the atomic displacement of the k-th atom in the Z-th 
unit cell given by: 


Uik = e k cos(q ■ Ri) q G {q L4 , q Mi } (5) 

where Rj is the Z-th lattice vector and e k gives the am¬ 
plitude of the displacement for the fc-th atom in the unit 
cell of the origin. 

If E(d) is the energy of the system per supercell as 
a function of the 2 x 2 x 2 distortion, in the harmonic 
approximation it is: 


E(d) ~ £7(0) + 1 £ 


d 2 E 


2 ddiaddjp 


dia djp 


d =0 


— £7(0) “ 1 “ ) a Cia,j0 dia djf3 


(6) 


(7) 


iaj/3 


and by grouping the two indeces (ia) = I we obtain a 
real-symmetric 72 x 72 matrix Cjj which has N = 72 
couples of real eigenvalues and eigenvectors (A, d^), so 
that for the distortion d ^ the system has the variation 
of energy: 


dE = ± ||d (A) || 2 ( 8 ) 

where || || is the euclidean norm of d^ x \ A negative 

eigenvalue corresponds to a displacement which lowers 
the energy of the system. 

By using the DFPT we calculated Cu and subse¬ 
quently we diagonalized it. Because of the symmetry, we 
obtained the same spectrum for the three spaces Vl* and 
the three spaces Vmi with the corresponding displace¬ 
ments related by threefold rotations. Consistently with 
the phonon analysis we found, for each of these space, 
two eigenspaces with symmetry A n , one of them with 
negative eigenvalue. Since we are interested in the insta¬ 
bilities of the system we focus on this kind of distortions. 


B. The A u distortions 


The displacements A u , for a point Li or M i: are 
transversal to the direction of propagation qr, planar (i.e. 
with zero component outside the xy- plane) and opposite 
for selenium atoms Sei and Se 2 on two adjacent wave- 
fronts. We indicate with —v the ratio between the dis¬ 
placements of the titanium atoms on a wave-front and 
the selenium atoms Sei on an adjacent wave-front: 


r e k ± q 

4 = 0 

CSe 2 = —€Sei 

, € T i = -z'esei 


(9) 


Thus, the sign of v indicates if the displacement of two 
adjacent Ti and Se front waves are in phase (y < 0) or 
out of phase (y > 0) (see Fig. |3|. 

These modes form a two dimensional vector space 2 A u 
made of displacements along a fixed line, where the only 
two free parameters left are the values of the shifts. We 
indicate with STi the shift of the Ti atoms on a dis¬ 
placement wave-front along a given direction and with 
£Se the displacement of the Se atoms on an adjacent 
wave-front along the opposite direction (see Fig. |3|. It is 
v = £Ti/£Se. By using the two parameters 5Ti, S Se the 
space 2 A u can be represented as in Fig.[4j By symmetry, 
two displacements with opposite values for both STi and 
£Se are equivalent whereas by changing the relative phase 
between the shifts of the Se and the Ti atoms we obtain 
two different configurations. In fact, a significant pa¬ 
rameter is the ratio v which identifies a one dimensional 
subspace A u {y) of 2 A u . Therefore, a general d G 2 A u is 
uniquely identified either by v and ||d|| or, for example, 
by v and £Ti (5Se). 

Each of the two orthogonal one-dimensional 
eigenspaces of symmetry A u found by diagonaliz¬ 
ing Cu is characterized by a specific value of the ratio, 
one positive and the other negative. In fact, it can 
be shown that two generic one-dimensional subspaces 
A u {y i) and A u { z/ 2 ), corresponding to different values of 
the ratio v = v\ and v = z/ 2 , are orthogonal if and only if 
^ 1^2 = — 2. We found that, in all the studied cases and 
in both the points M and L, the A u with the smallest 
eigenvalue A corresponds always to the positive ratio 
v > 0, that is to the out of phase distortion. The values 
found for the smallest A and the corresponding v in the 
point L (Al, vl) and M (Am, Fm), for all the cases, are 
reported in Tab. El From now on we analyze only the 
cases where the system displays instability, so we do not 
consider LDAxh anymore. 

The value vl (fm) defines (up to a sign) three orthog¬ 
onal unit vectors (c^mJ corresponding to degenerate 
displacements of type A u which generate a three dimen¬ 
sional space (V^). By considering displacement vec¬ 
tors d in this space with an increasing moduls ||d|| but 
fixed direction, we can study the energy variation of the 
system along a pattern. We observe that the energy, after 
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9 Ti • Se 

V > 0 

- A u (XiandMi) 



v < 0 



A u (I 3 andM 3 ) 



3gjr and 3 q M 



FIG. 3. (color online) Schematic representation of the atomic displacements in the 2x2x2 cell of a single layer viewed from 
above for the three A u modes in the points Li and Mi and the triple modes 3 q L and 3 q M . Both the not equivalent cases for 
v > 0 (left-hand figures) and v < 0 (right-hand figures) are shown, v being the ratio between (5Ti and #Se. The Se atoms with 
a dot are on the upper plane with respect to the Ti atoms plane. In the first picture it is also drawn the direct lattice basis a* 
and the lxlxl unit cell. The two central arrows indicate the positive directions used to measure the shifts £Ti and $Se. In 
the 3 q mode with i/>0a three atom cluster is highlighted with a dotted line. Notice that the triple-g modes shown are not 
equal to the sum of the A u single-g modes represented as, otherwise, the displacements should be larger (the displacements of 
the atoms in a triple mode are two times larger than the displacements of the component single modes). 
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2 A u 



V2 6Se (A) 


FIG. 4. (color online) Diagram representing the distortion 
space 2 A u in L and M as a function of the displacements 
(£Se, 5Ti). The factor y/2 on the horizontal axis comes from 
the two Se atoms in the unit cell and it is necessary in order 
to convert the orthogonality condition in 2 A u into the eu¬ 
clidean orthogonality on the diagram. A line on this plot rep¬ 
resents a one dimensional subspace A u (y) with a specific ratio 
v — STi/SSe. As an example, the two orthogonal vectors cor¬ 
responding to the A u distortion modes of the GGA^h W case 
are drawn. The point corresponding to the CDW distortion 
experimentally observed^ is also drawn. 


TABLE IV. Values of the eigenvalue A and the ratio v — 
5Ti/$Se for the A u distortions in L and M corresponding to 
the smallest eigenvalue. Last row: value of the ratio measured 
for the CDW distortiorP. 



Am 

vm 

Al 

vl 

LDAexp 

-1.07 

2.89 

-1.37 

2.93 

GGAexp 

-0.89 

2.60 

-1.26 

2.53 

ggaLT 

-0.84 

2.64 

-1.20 

2.56 

LDAxh 

+0.69 

2.34 

+0.25 

2.42 

GGAxh 

-2.26 

2.46 

-2.25 

2.43 

GGAx 3W 

-1.14 

2.51 

-1.31 

2.46 

Exp 

Tdw 

- 

- 

- 

3.42 + 1.48 


an initial parabolic decrease, starts departing from the 
harmonic regime, reaches a minimum and then increases 
(c/. Fig. |J. The minimum along this energy path cor¬ 
responds to the configuration which gives the most sta¬ 
ble structure obtainable with that kind of distortion. In 
general, different patterns in (V^) return different 
results as for finite displacements the symmetries of the 
lattice are not preserved. 


C. The single-point and triple-point patterns 


The displacement pattern of type A u characterized by 
a unit vector d/,. (cZmJ is also called a single-q L pattern 
( single-q M pattern). In Fig. [5] we show, for the several 
studied cases, the energy path for the single-^ patterns 
in L and M. As we can see the distortion of type L re¬ 
turns always a structure more stable than the distortion 
of type M, except in the GGA T h case where the two en¬ 
ergy patterns approximatively coincide, meaning that the 
interaction between layers plays a not negligible role in 
the minimization of the total energy unless their distance 
is large enough. Moreover, as long as the cell dimensions 
are comparable, we find similar results irrespective of the 
local functional used. Instead in the GGAxh case, where 


the volume of the unit cell is larger (see Tab. VI), the 
energy gain due to the distortion is higher. These con¬ 
siderations leads to the conclusion that the suppression 
of the instability in the LDAxh case is ascribable to a 
pure volume effect, since the LDA theoretica l un it cell is 
smaller than the experimental one (see Tab. VI). 

The pattern characterized by the unit vectors d^L and 
d 3 M obtained by combining the three directions d/,. and 

d Mi 


3 1 

dsM = 

V 3 



( 10 ) 


are called a triple-q L (3 q L ) pattern and a triple-q M 
(3<Zm) pattern, respectively. In fact, by definition, a gen¬ 
eral triple-^ distortion of type L (M) is obtained by su¬ 
perimposing, with equal weights, three A u distortions for 
the points Li (Mi) having the same values of £Ti and 5Se. 

In a layer, the shifts of the atoms in a 3 q distortion 
can be described by considering the TiSe6 octahedral 
structure of the systeirPS(see Fig. [I]). We distinguish two 
kinds of Ti and Se atoms: in one the Ti(a) atoms do not 
move and are in the middle of an octahedron Ti(a)Se(a)6 
where the three Se(a) atoms above Ti(a) and the three 
Se(a) atoms below Ti(a) rotate with opposite direction. 
Thus, there are couples of Se(<a)’s which become closer 
and, depending on weather the component three modes 
have v > 0 or v < 0, they attract or repulse a Ti (/3) atom 
(in the first case, as a consequence, we observe a Ti(/?)- 
Se(a) bond shortening and the formation of three-atoms 
clusters Ti(/J)Se(a )2 in the system). In the distortion we 
also have S e(/3) atoms which are not involved in any ro¬ 
tation (they are originally in octahedra centered around 
Ti (/3) atoms) but stay in their position. 

For an adjacent layer, depending on whether we are 
considering a 3 q M or a 3 q L distortion, the displacement 
of the atoms is the same or the opposite one ( i.e . in the 
3 q L case if in one layer the Se(a)’s on the upper plane 
and lower plane show a clockwise and a counterclockwise 
rotation the opposite happens in an adjacent layers, re¬ 
spectively). Two 3 q distortions with v > 0 and v < 0, in 
a layer and view from above, are shown in Fig. [3j 
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Mil (A) 


FIG. 5. (color online) Variation of the energy with respect to 
the undistorted phase obtained by moving the atoms accord¬ 
ing to different patterns. The unit vectors 
characterize a single q L and q M . mode, and a triple q L 
and q M mode, respectively. The vertical dashed line in the 
GGA^h W plot marks the point where the calculations of Fig. [7] 
have been performed. 


D. The CDW distortion 

In their seminal paper Di Salvo and coworkers stated 
that, according to neutron-diffraction measurements and 
symmetry considerations, the lattice distortion experi¬ 
mentally observed with the CDW is a 3 q L mode with 
v > 0, with data taken at 77 K which best fitted with 
the valued 1 : 


£ Ti Ex P = (0.042 =b 0.007)A 

<5Se Exp = (0.014 ± 0.004)A '' 11} 

for the displacements of the atoms in the single Li com¬ 
ponent modes. These values correspond, for the complete 
3 L pattern, to the displacements: 


£ Ti Ex P = ( 0>085 ± 0> oi4)A 

(3) 

£Se Exp = (0.028 ± 0.007)A 


( 12 ) 


for the Ti and Se atoms which actually move (in a 3 q 
mode not all the atoms move). The experimental values 
for the displacements in a single q L component mode 
are also shown on the diagram in Fig. [4] and correspond 
to the experimental estimate for the ratio: 

z/ Exp = <m Exp /£Se Exp = 3.42 ± 1.48 (13) 


Motivated by these experimental results we calculated 
the energy pattern of the 3 q L mode in for the sev¬ 
eral studied cases. Moreover, in order to study the role 
played by the interaction between layers, we also calcu¬ 
lated the energy pattern of the 3 q M distortion in 
The results are shown in Figj5] As we can see a triple -q 
pattern returns always a structure more stable than the 
corresponding single-^ displacement and it is always the 
3 q L distortion wich gives the lowest energy (except in 
the GGA T h case where the modes 3 q L and 3 q M are al¬ 
most degenerate). Moreover, as for the single-q patterns, 
when the cell dimensions are comparable we find similar 
results, irrespective of the local functional used. 

In Tab. [V] we report, for the minimum point of the 
3 q L energy pattern, the values of the energy variation 
(per supercell) and the shift of the atoms (for the com¬ 
ponent A u single-qr L modes) with respect to the undis¬ 
torted crystal. The values of the displacements for the 
several cases, compared with the experimental result, are 
also shown in Fig. [ 6 ]( which corresponds to a portion of 
the diagram in FigTHl). 

The fact that a 3 q L pattern gives a structure more sta¬ 
ble than the single q L pattern is in agreement with the ex¬ 
perimental findings but, in order to demonstrate that the 
ab-initio calculations are able to predict the experimental 
results, it is necessary to go further and demonstrate that 
the 3 q L pattern gives the most stable structure among 
all the possible distortions in V£. Since we expected 
analogous results in all the cases analyzed we just con¬ 
sidered the GGAxh W one. In order to accomplish the 




































TABLE V. First column: largest energy gain (per supercell), 
with respect to the undistorted phase, for the 3 q L triple pat¬ 
tern in V£\ Second and third columns: corresponding atomic 
displacements, with respect to the undistorted phase, for the 
component single modes (A u symmetry). Notice that for the 
resultant 3 q L mode the displacement of the atoms which ac¬ 
tually move is two times larger, cf. Fig. [3] First row: experi¬ 
mental measure of the displacement for the CDW phase with 
respect to the high-temperature phasJ^. 



AF min (meV) 

<m(A) 

<5Se(A) 

EXP 

- 

0.042 ±0.007 

0.014 ±0.004 

LDAexp 

-8.4 

0.030 

0.010 

gga Bxp 

-7.4 

0.027 

0.011 

GGAgw 

-6.8 

0.026 

0.010 

GGAth 

-23.4 

0.037 

0.015 

GGAxh W 

-8.3 

0.028 

0.011 


the sphere of radius V to obtain a general scan of the 
surface: 

d = t c,d L , f> |!=P (14) 

i' = 1 r\ 

Ci> 0 

An heatmap plot of the results obtained is showed 
in Fig.[7| our calculations confirm that on this sphere the 
3 q L pattern (which has components q = V/y/3) returns, 
among all the possible patterns in , the most stable 
structure. We conclude that, in the frame of the electron- 
phonon interaction, by using first-principle calculations 
we are able to recover the CDW structural instability 
experimentally observed for TiSe 2 - 


2 A u 
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0.04 


0.03 




0.02 


0.01 



0.01 0.02 

V2 6Se (A) 
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FIG. 6. (color online) Filled points: displacement of the 
atoms corresponding to the energy minimum along the 3 L 
pattern in V f (single A u component). Empty points: A u 
component of the total displacement obtained by further re¬ 
laxing the structure. 


task, we should have considered a uniform sample of the 
unit sphere in Vf and computed the energy path for an 
increasing modulus of the distortion along each of these 
directions. Instead, in order to reduce the workload, we 
just computed the energy of the system for distortions d 
having a fixed modulus ||d|| =T>, where V = 0.146A is a 
length approximately in half position between the mini¬ 
mum along the q L . and the 3 q L patterns (cf. Fig. |5]). In 
this way we only needed to scan the energy of the system 
for a uniform grid on the two dimensional sphere in Vf 
with radius V. We used the basis di Ji to parametrize 
the space Vf and the fact that given a general vector 
d = i c i&Li all the other vectors obtained from it 
by changing the sign of its components q give equivalent 
distortions. Therefore, we considered only one octant of 


E. Optimization of crystal structure in the 
CDW phase 


In order to find the equilibrium configuration for the 
distorted system we subsequently relaxed the structure 
starting from the configuration corresponding to mini¬ 
mum of the energy along the 3 q L path. We relaxed the 
whole structure (cell and internal positions) or only the 
internal positions depending on whether we were con¬ 
sidering the theoretical or the (fixed) experimental cell. 
This led to a further small gain in en ergy with respect to 
the undistorted phase (see Tab . |vT|) . 

As it can be seen from Tab. |VII[ the effect of the re¬ 
laxation on the lattice parameters is quite small, with 
a relative variation of the order of 10 -3 . For the in¬ 
ternal displacements we define the vectors in V, d mm 
and d rlx , characterizing the atomic shifts from the undis¬ 
torted phase to the minimum along the 3 q L path in Vf 
and the equilibrium position reached after the relaxation, 
respectively (in the theoretical cell cases d rlx is the rel¬ 
ative displacement with respect to the cell). A quan¬ 
titative measure of the extent of the additional atomic 
displacements due to the relaxation is obtained by com¬ 
paring the modulus and the direction of these vectors. 
The results are shown in Tab. m and, as we can see, 
the two vectors are very similar (only slightly different in 
the GGATh case), meaning that the relaxation does not 
give an additional huge displacement of the atoms. 

The displacement of the atoms due to the optimization 
is obviously restricted by the symmetry of the distorted 
phase. The total displacements d rlx is mainly made of a 
3 q L distortion but it also has a small component which 
changes the value of R shown in Fig. [l] (the value of this 
parameter being not fixed by the symmetry, anymore) 
and modifies the distances h a and hp between the not 
equivalent Se(a) and S e(/3) atoms and the Ti plane. As 
a consequence, in the distorted structure the upper (and 
lower) Se atoms of a layer are not on the same plane 
anymore (cf. also Ref. f26[) . The values found for the 


atomic displacements are reported in Tab. VI and, in 
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FIG. 7. (color online) GGA^h W case: heatmap plot of the energy variation with respect to the undistorted phase for distortions 
corresponding to vectors in V^ having modulus V — 0.146 A. Due to the symmetry of the system the result is shown on a 
single octant. Left hand panel: 3D color plot on the spherical surface of radius V. Right hand panel: 2D plot in spherical 
coordinates ($, c p) associated to the basis (gRJ. The lowest energy is obtained for the 3 q L distortion, whose angular coordinates 
are Lp — 7r/4, $ = arccos(l/v / 3)- 


TABLE VI. Effect of the relaxation on the minimum of the energy along the 3 q L path in Vtotal variation of the energy (per 
supercell) from the undistorted phase AE rlx , shifts £Ti and (5Se for the A u component of the atomic displacements, variation 
of R (see Fig. [l]), and change of h for the two not equivalent atoms Se(a) and Se(/3). First row: experimental measure of the 
displacement for the CDW phase with respect to the high-temperature phasj^. 



A£ rlx (meV) 

<5Ti(A) 

5Se(A) 

SR/R 

5h a /h 

5hp/h 

EXP 

- 

0.042 ± 0.007 

0.014 ±0.004 

~ 

- 

- 

LDAexp 

-9.4 

0.030 

0.010 

-0.0003 

0.0011 

-0.0001 

gga Bxp 

-8.7 

0.028 

0.011 

-0.0003 

0.0016 

0.0004 

GGA^w 

-7.9 

0.027 

0.010 

-0.0002 

0.0014 

0.0004 

GGAxh 

-27.6 

0.040 

0.016 

-0.0018 

0.0031 

-0.0038 

GGA™ W 

-9.5 

0.029 

0.011 

-0.0004 

0.0010 

-0.0006 


TABLE VII. Effect of the relaxation on the minimum of the 
energy along the 3 q L path in : variation of the cell param¬ 
eters a and c (in the theoretical cell cases) and comparison of 
the vectors d mm and d rlx representing the shifts of the internal 
atomic positions from the undistorted phase to, the minimum 
of the energy along the 3 q L path and the final relaxed con¬ 
figuration, respectively. 



l|d rlx || 

(d rlx d min ) 

A a 

Ac 


||d mi ”|| 

||d rlx || ||d™ || 

a 

c 

LDAexp 

1.0008 

0.9992 

- 

- 

GGAexp 

1.0016 

0.9981 

~ 

- 

GGAgw 

1.0022 

0.9985 

~ 

- 

GGA T h 

1.0392 

0.9918 

0.0006 

0.0013 

GGAxh W 

1.0046 

0.9990 

0.0006 

0.0006 


In conclusion, the energy gain and the A u component 
of the atomic displacements obtained in all the studied 
cases are similar for similar unit cell and the atomic shifts 
are not far from the experimental measures. In par¬ 
ticular, the ab-initio calculations based on the electron- 
phonon interaction predict the correct displacement pat¬ 
tern for the low energy distortion with the values of the 
atom displacements slightly underestimated (especially 
for the Ti). Only for the GGAxh case, which however 
uses cell parameters quite different from the experimen¬ 
tal ones, we have found different results with, in partic¬ 
ular, a larger atomic displacement and a larger energy 
decrease for the low temperature phase with respect to 
the other cases. 


particular, the updated values of c)Ti and £Se for the A u 
component of the displacement are also shown in Fig. [6] 
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F. Chiral Charge-Density Waves 

Quite recently, on the basis of tunneling microscopy 27 
and x-ray scattering®! measurements, it has been re¬ 
ported the existence of a second phase transition for 
TiSe 2 which follows the conventional CDW transition 
and occurs at temperature around 7 K below the Tcdw- 
As we have seen, the conventional CDW transition is due 
to the superposition with equal weights of three modes 
describing a charge transfer process from Ti-3d to Se-4p 
orbitals and a related commensurate lattice distortion. 
The new type of transition is still given by the sum of 
these three charge transfer waves but with different rel¬ 
ative phases, producing an helical distribution of charge 
along the z axis of the 2x2x2 cell. As a consequence, 
there are two possible, energetically degenerate, distribu¬ 
tions of charge characterizing two chiral phases related by 
a mirror transformation. These are the so called Chiral 
Charge-Density Waves (CCDW) and a driving mecha¬ 
nism for chiral symmetry breaking it has also been pro¬ 
posed. 29 One of the most important features of CCDW 
transition is the reduction of the crystal symmetry from 
a threefold to a twofold rotation symmetry. 

In our simulations we did not see any evidence for a 
phase different from the conventional CDW: the lowest 
total energy for the system was obtained by consider¬ 
ing the superimposition of the three modes with equal 
weights. Nonetheless, we find worthwhile to stress that, 
by using a very simple argument, it is still possible to ex¬ 
plain the existence of different chiral domains on the sur¬ 
face of a TiSe 2 sample (but not the change of symmetry). 
As shown in Fig. [TJ the CDW transition is characterized, 
in a cell, by a clockwise or counterclockwise rotation of 
the Se atoms in the upper layer (and the opposite for the 
Se in the bottom layer). Thus, for a cell, there are two 
possible degenerate but different CDW structures that 
cannot be superimposed solely with rotational transfor¬ 
mations and this could naturally led to the formation, 
on the surface of a TiSe 2 sample, of two different chiral 
domains. We think that this aspect should be properly 
taken into account when the results of an experiment 
about the CCDW is interpreted. For the same reason, 
we believe that a simple surface experiment, like a scan¬ 
ning tunnel microscopy measure, could be misleading and 
that only a measure of the symmetry change of the crys¬ 
tal can prove the existence of a second phase transition 
different from the conventional CDW. 


IV. ELECTRONIC STRUCTURE IN THE 
CDW PHASE 

A. Undistorted energy bands 

In this section we discuss the evolution of the DFT 
band structure under the lowest energy distortion found 
in V^. In Fig. [8] we can see the undistorted bands of 
TiSe 2 around the Fermi level (Ep) for an high-symmetry 


path (c/. Ref. [30]). The atomic/orbital character of the 
bands is also shown (c/. Ref. [26]). In all the cases we find 
similar results: around the Fermi level there are only Ti- 
3 d and Se-4p derived bands with a narrow Ti-3d band 
which is almost entirely unoccupied except around the L 
point where it crosses the Fermi level and forms a hole 
pocket; moreover, two Se-p and Ti -d strongly hybridized 
bands cross the Fermi level. A qualitative difference is 
instead found for a narrow Se-4p band whose position 
mostly depends on the cell dimensions (in particular on 
the distance between the Se and the Ti atom planes) but 
not on the local functional used. In fact, as we increase 
the value of c we observe a lowering of the band max¬ 
imum in T (in addition to an increase of the structural 
instability, as we have seen in the previous section). Nev¬ 
ertheless, as we will see, this band does not change during 
the distortion. 


B. Energy bands and DOS in the CDW phase 

We are interested in the evolution of the band struc¬ 
ture induced by the 2x2x2 distortion. From now on 
we consider the cases LDAe xp and GGA^ W . With the 
subscript ‘sc’, set below the labels T, M, K of the special 
points in the BZ, we refer to the corresponding points of 
the 2 x 2 x 2 super-cell Brillouin zone (SBZ) which has 
the same shape of the BZ but half size. 

In Fig.[9]we show the 2x2x2 bands, around the Fermi 
level, for an high symmetry line of the SBZ for the undis¬ 
torted and distorted phase, the last one corresponding to 
the minimum of the energy along the 3 q L distortion path 
in . On formation of the superstructure we observe a 
neat change near the Fermi level with similar character¬ 
istics in the two cases. In particular, around T sc we find 
avoided crossings at the Fermi level and, just below it, the 
appearance of a (reversed) Mexican-hat structure due to 
the repulsion of bands having Ti-3d character (see Fig. [9]) . 
The Se-4p bands crossing the Fermi level around the wave 
vector (1/2)T SC M SC , instead, do not suffer any modifica¬ 
tion during the distortion, as anticipated. As expected, 
the change of the electronic dispersion around the Fermi 
leads to a change of the Density of States (DOS) as it is 
shown in Fig. [lO] On forming the superlattice structure 
the DOS decreases by around 40% at Ep and a peak, 
essentially due to the Ti-4d orbitals, develops at 0.15/0.2 
eV below Ep. These effects have been qualitatively ob¬ 
served in some previous theoretical works (see Ref. [31] 
for example, for a tight-binding study) and are, at some 
level, compatible with the change of the resistivity ex¬ 
perimentally observed in the CDW transition 1 . 


C. Energy bands folding and unfolding 

The superlattice distortion doubles the original lattice 
periodicity of the system. In fact, each eigenfunction 
^k,j of the distorted system has the pseudo-momentum 
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FIG. 8 . (color online) Band structure of TiSe 2 (undistorted phase) for several cases along a BZ high-symmetry line ( cf . Fig. [ 2 ]). 
The atomic/orbital character of bands is expressed by using different colors. Starting from the upper left-hand figure, the 
dimension of the cell along c increases clockwise. Only in the LDATh case the system does not have an instability. 


K in the SBZ, which is one-eighth of the original BZ. 
Nevertheless, the function ^k,j is made of eight con¬ 
tributions V’fc.’ which are functions pseudo-periodic on 
the original lattice and whose pseudo-momenta ki are 
obtained by unfolding K to the original BZ: 

^k,j = X! <’ J ( 15 ) 

i=l 

The spectral weights c^ ,J : 

<’ j «ii<’ j ii 2 E<’ J = i as) 

i=1 

can be used to evaluate the contributions to ^k,j coming 
from different points of the original BZ (see appendix |A|). 


From geometrical considerations we see that in our case 
the SBZ can be unfolded into eight regions of the BZ cen¬ 
tered around the points T, A, Li, Mi, respectively, and we 
can use this property to label the corresponding unfold¬ 
ing weights. Moreover, due to the threefold symmetry 
of the system, for a point K it is convenient to sum the 
contributions coming from the three equivalent Li points, 
and the same for the contributions coming from the three 
Mi points: in this way we have four contributions of type 
T,L, M, A depending on the BZ portion they came from. 
In Fig. [TT] the weights of these four contributions on the 
bands are shown by means of a color code. In this way 
we easily recognize, for example, that the characteristic 
band configuration under Ep in r sc has a pure L char¬ 
acter. 

A complementary method to describe the 2 x 2 x 2 
distortion in the frame of the original translation sym- 
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FIG. 9. (color online) DFT bands, with orbital character, of the 2x2x2 super-lattice along the high symmetry line r s cM S ciF sc r s c 
of the super-cell Brillouin zone. Left-hand panels: undistorted configuration (||d|| = 0.00 A). Right-hand panels: distorted 
configuration (||d|| = 0.16 A). An arrow highlights, in two cases, the effect of the repulsion between two Ti-3d bands. 
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E-E f (eV) 


FIG. 10. (color online) Density of states around the Fermi 
level for three cases in the undistorted (red line) and distorted 
(green line) phase. 


metry (and best suited to make a direct comparison with 
ARPES experiments) is to unfold the superlattice band 
structure from SBZ to the original BZ 32 33 . This consists 
in plotting for the points k of a line in the BZ the energy 
bands Ej^ j of the distorted system with an intensity 
J ec l ua l to the spectral weights (for = 0 

we have full transparency, i.e. no band, and for J = 1 


full opacity). In this way the displacement of the atoms 
is observed as a distortion, smearing and fade of the orig¬ 
inal bands, plus the appearance of new ghost bands. In 
Fig. [12] the unfolded bands with the orbital character are 
shown. In order to ease the comparison with the undis¬ 
torted phase we also superimpose the unfolded bands on 
the original band structure. From these figures we can 
clearly see what is principal effect of the 3 L distortion: 
the Ti-4d bands which were crossing the Fermi level now 
open a gap whereas the Se-3p bands essentially remain 
unaffected. Moreover we see the appearance of a ghost 
band in L which gives the characteristic structure we 
have already discussed. 


V. COMPARISON WITH ARPES 
A. LDAexp and GGA™ 

In the previous sections we analyzed the kind of distor¬ 
tion which lowers the total energy of the system and the 
consequent changes in the electronic structure as they 
are found by ab-initio DFT calculations. In this section 
we compare the results of the calculated band structure 
and the ARPES data taken from Ref. |34j In Fig. [13] we 
show the DFT bands for the undistorted and distorted 
phases superimposed on the ARPES data taken at high 
(T = 300 K) and low (T = 35 K) temperature, respec¬ 
tively. Due to the indetermination of k z in the ARPES, 
and the substantial k z dispersion in the electronic struc¬ 
ture, we consider the DFT bands along the two directions 
T — M and A — L in BZ. In the high temperature case we 
simply plot the bands along these two lines for the undis¬ 
torted system and in the low temperature case we plot 
the 2 x 2 x 2 bands of the distorted structure unfolded 
to BZ along these two directions. In this second case, in 
order to ease the comparison with the underlying figure, 
we slightly enhanced the intensity of the unfolded bands 
by scaling the whole transparency by a factor / = 4: 

n = (17) 

As we can see, the effect of the distortion on the the¬ 
oretical bands can be considered in reasonable agree¬ 
ment with the experimental data for the low temperature 
phase (in particular it is remarkable that the Mexican- 
hat found in L seems to reproduce an experimentally 
observed feature). Nevertheless, there are features in the 
bands which seem to be not compatible with the ARPES 
findings. This is even more evident for the high temper¬ 
ature phase since in this case the ARPES data are not in 
good agreement with the bands of the undistorted struc¬ 
ture. In the second case, a possible reason for the mis¬ 
match could be the effect of phonon quantum fluctuations 
on the electronic structure, which are obviously stronger 
when we consider higher temperatures and which are not 
taken into account in our “zero temperature” calcula¬ 
tions. More generally, a possible explanation for the not 
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LDA Exp + *7 = 3.9 



FIG. 11. (color online) DFT bands of the 2x2x2 super-lattice along the high symmetry line r sc M sc lF S cr s c of the super¬ 
cell Brillouin zone. Left-hand panels: undistorted configuration (||d|| = 0.00 A). Right-hand panels: distorted configuration 
(||d|| = 0.16 A). The colors indicate the weights of the corresponding eigenfunctions on different parts of the lxlxl Brillouin 
zone. 
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FIG. 12. (color online) Unfolded bands for three cases. Left-hand panels: unfolded bands for the distorted phase on a high 
symmetry line of the original Brillouin zone with the orbital characters highlighted by colors. Right-hand panels: unfolded 
bands (in red) superimposed on the undistorted bands (in black) along the same line. 
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perfect agreement between ARPES and the calculated 
bands in both the phases could be the effect of the cor¬ 
relation due to the localized d orbitals of Ti (correlation 
which is not taken into account in a pure DFT calcula¬ 
tion). We explored this possibility by using the DFT+U 
method (in particular the LDA+J7 method) and in the 
next section we show the results of this analysis. 


B. LDAexp+H 


As explained in Sec. [TT| we considered the Hubbard¬ 
like correction to the electronic structure of LDAexp- We 
used the experimental cell, the internal theoretical coor¬ 
dinates obtained with LDA and, on the top of this, the U 
correction for the Ti -d orbitals. The most evident effect 
of introducing U is to open a gap between the bands with 
the result of obtaining, for U ~ 3.8 eV, a metal-insulator 
transition (see Fig. 14). Motivated by this result, we 
found appropriate to estimate ab-initio the proper value 
of U by using a self-consistent procedure. In a few steps 
we obtained the converged value U = 3.902 eV for which 
the system is in an insulating phase with a small gap 
of approximatively 0.014 eV (see Fig . 14). This result 
is in line with a recent observation^, although another 
optical experiment 36 ) has found a semi-metallic state in 
both the high-temperature and low-temperature phases. 
Nevertheless, the LDAexp +U result seems to be closer 
to the experimental observations than the large negative 
gap obtained with pure DFT calculations. 

We used the same value of U also for the distorted 
structure found in the LDAexp case. We observe that the 
distortion has, on the electronic bands, the same quali¬ 
tative effect already observed without 1/, the shift of the 
Ti -d bands now increasing the initial small band gap up 
to a value of approximatively 0.2 eV (see Fig. Emm 


12). In the lowest panel of Fig. [I3lwe compare the cal¬ 


culated LDAexp + U bands with ARPES data: now the 
theoretical results for both the distorted and undistorted 
structure are in very good agreement with experimental 
data taken at high and low temperature, respectively. 

According to these results, the correction provided by 
U seems to solve even the last dubious results regard¬ 
ing the comparison between the DFT electronic struc¬ 
ture and the experiment, in particular by giving a per¬ 
fect match between the theoretical bands and the ARPES 
data. Unfortunately, a serious drawback of this approach 
is that the presence of U eliminates the instability. 

In Fig [T5] we show the lowest phonon frequency ujl 
(ujm) in L (M) as a function of U calculated by us¬ 
ing the finite-difference method. When a phonon fre¬ 
quency uo is imaginary we conventionally indicate it with 
the negative value — |cj|. We considered the LDA system 
with experimental cell and internal coordinates now ob¬ 
tained by relaxing the atomic positions for each value of 
U. It can be observed that the frequencies, which are 
negative (i.e. imaginary) for U = OeV, increase as we 
increase U and begin to have a positive value around 


U ~ 2.5 eV for ujl and U ~ 1.5 eV for ujm- A phonons 
converged calculation with this method was revealed to 
be extremely costly, especially for values of U close to 
the metal-insulator transition, so we also explicitly veri¬ 
fied that the total energy of the system increases as we 
move the atoms along the expected distortion pattern. 
We performed similar tests also for the other cases and 
we reached a similar conclusion: the DFT+U system is 
not unstable anymore for the value of U which seems 
to return the best agreement between the undistorted 
bands and the high temperature ARPES data (notice 
that this value is functional-dependent, U being not a 
physical measurable quantity). 

At this point an intriguing question arises: on one hand 
it seems that we are able to correct the energy band 
structure by using [/, but on the other hand the same 
correction spoils the prediction of the instability driven 
by the electron-phonon coupling, a prediction whose re¬ 
liability, in turn, strongly depends on the correctness of 
the energy evaluation. 

A first possible explanation for this apparent paradox 
could be in the simplistic approach used to describe the 
structural instability with U. In general, the energy of a 
DFT+U system (with cell fixed) is given by a function 
(o(Ri,Uj) of the atomic positions Ri and the values Uj 
for the orbitals. Moreover, for a certain atomic configu¬ 
ration Ri we have a proper value Uj = %j(Rj) for the 
Hubbard terms. So the total energy of the system as a 
function of the atomic positions, E(Ri ), is given by: 


E(R i )=g(R i ,%(R i )) (18) 

and when we consider the forces and the force constant 
matrix (through the first and second derivative of E with 
respect to the positions R^ respectively) we should take 
into account the full dependency on the atomic positions. 
On the contrary, a commonly accepted procedure is to ne¬ 
glect the dependence of the Hubbard terms on Ri and to 
move the atoms by keeping fixed the values of Uj found 
for the undistorted system. This is the approach we con¬ 
sidered in our study but its correctness seems to be ques¬ 
tionable in this case since even a small variation of U 
(so small that it does not affect the energy bands po¬ 
sition) leads to a total energy variation which is huge 
if compared with the energy decrease obtained with the 
distortion. For example, during the self-consistent calcu¬ 
lation of U, we found that for AU ~ 1 meV the system 
has a variation of the total energy A E ~ 8 meV. 

Therefore, in general, we expect that in order to make 
a proper structural analysis of TiSe 2 in the DFT+U 
scheme it is necessary to take into account the full posi¬ 
tion dependence (and consider that in the distorted phase 
we have two not equivalent Ti sites with, in principle, 
two different values of U). But it also means that, for 
example, a careful study of the convergence of U with 
the size of the supercell used to perform the perturbative 
calculation should be done. Nonetheless, it seems that, 
in general, the precision required for the value of U (at 
least 10 -5 eV) is unrealistic and probably rules out the 
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FIG. 13. (color online) Comparison with ARPES data (from Ref. f34l) in three cases. First row: ARPES data taken at high 
temperature (left-hand panel, T = 300 K) and low temperature (right-hand panel, T = 35 K). Last three rows: DFT bands 
superimposed on the ARPES data. Black lines: F — M direction in BZ. Red lines: A — L direction in BZ. In the high temperature 
cases the bands of the undistorted structure are plotted. In the low temperature cases the 2x2x2 bands of the distorted 
structure unfolded to the original BZ are plotted. The intensity of the unfolded ghost bands have been slightly enhanced in 
order to ease the comparison with the ARPES figure (see main text). 
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U (eV) 


FIG. 14. (color online) Effect of U on the bands. In the last 
figure the vertical dashed line indicates the final self-consistent 
value U ~ 3.9 eV, which corresponds to an insulating phase 
with a small gap of approximatively 0.014 eV. 


use of a DFT +U method in this form. 

A second consideration relies on the variation shown 
by the frequencies ujl{U) and ujm{U) with U. In fact, 
we observe a monotonic increase of these two quantities 
until the points A p and Y p in Fig. 14 go below the point 
L p and the point M p , respectively; after that the value of 
the frequencies essentially remain constant. So it seems 



0.0 1.0 2.0 3.0 4.0 5.0 

U (eV) 


FIG. 15. (color online) Lowest phonon frequency in L and M 
as a function of U, for the LDA -\-U system with experimental 
cell and theoretical internal coordinates. 


that U removes the instability because of the change in 
the relative positions of the bands. Thus, it could be 
impossible to maintain the lattice instability and, at the 
same time, reproduce the high temperature ARPES ex¬ 
periment. This would confirm the difficulty of predicting 
the high temperature energy dispersion with a DFT(+C/) 
band calculation. 


VI. CONCLUSIONS 

In this work we have presented a first-principle study of 
the charge-density wave in TiSe 2 by using first-principle 
DFT calculations. We have considered several local func¬ 
tionals and both experimental and theoretical cell param¬ 
eters. 

We have shown that the results are robust with re¬ 
spect to the local functional used, whereas they depend 
essentially on the distance between the layers, a larger 
distance corresponding to an higher instability. In fact, 
no instability is observed in the LDA case with zero the¬ 
oretical pressure, in this case the distance between the 
layers being smaller than the experimental one. On the 
contrary, as long as the cell parameters are in agreement 
with the experiments, or even larger, the calculations re¬ 
produce not only the structural instability but also the 
triple-qr L distortion observed in neutron diffraction ex- 
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periments. We have also explicitly analyzed the role of 
the interaction between layers by considering distortions 
of type M and L. We have concluded that, in both the 
points of the BZ, is the triple-qr distortion which always 
gives the most energetically favorable configuration, elu¬ 
cidating the importance of the in-plane displacements in 
order to lower the total energy. Nevertheless, as long as 
the layers are close enough to have a significant interac¬ 
tion, the L distortions are always more favorable than the 
M ones, meaning that the intra-layer interaction plays a 
not negligible role in the CDW transition. 

We have also analyzed the changes in the electronic 
structure of the system with the CDW formation. In 
the experiment, the phase transition is accompanied by 
changes in the transport properties. We have found mod¬ 
ifications in the calculated bands near the Fermi level, 
with a resulting depletion in the DOS, which are com¬ 
patible with the changes measured in the resistivity. The 
bands for the superstructure have also been unfolded 
into the undistorted BZ in order to analyze their evo¬ 
lution with the phase transition and compare them with 
the findings of an ARPES experiment. We have shown 
that the unfolded bands of the distorted structure are in 
reasonable agreement with the low-temperature ARPES 
data. Nevertheless, the bands for the undistorted phase 
are in not good agreement with ARPES. 

In order to correct the mismatch between ARPES 
data and band calculations we have explored the role 
of the correlation of the Ti -d electrons by performing 
an LDA -\-U calculation. We have estimated ab-initio, 
self-consistently, the value of the parameter U for the 
localized Ti -d orbital in the LDAe xp undistorted struc¬ 
ture. For this value of f/, the Hubbard-like correction 
opens a small indirect gap in the electronic bands of the 
undistorted phase, whereas a pure DFT calculations pro¬ 
vides a semi-metallic configuration with a large negative 
gap. This result is compatible with some recent experi¬ 
ments. More importantly, the LDA +U bands are in very 
good agreement with the results of an high-temperature 
ARPES experiment. By using the same value of U for 
the distorted structure, we have shown that the CDW 
increases the size of the band gap. The unfolded bands 
of the distorted phase are also in good agreement with 
the low-temperature ARPES data. 

The drawback of the LDA -\-U approach is that the U 
correction removes the instability, the phonon frequency 
in M and L becoming real as we increase the value of U . 
A possible explanation for this effect has been identified 
in the neglected, but in principle existing, dependence of 
the value of U from the atomic positions. On the other 
hand, the loss of the instability could be the correct result 
of the electronic bands displacement caused by the local- 
chemistry correction given by U in TiSe 2 - This would be 
compatible with the suppression of the CDW observed in 
the ternary crystals TiS a ,Se 2 - ;E for x > 0.95 1 . In fact, at 
variance with TiSe 2 , TiS 2 has no CDW instability, but 
the electronic structures of these two systems are very 
similar, the relative position of the valence and conduc¬ 


tion band being the main relevant difference: while TiSe 2 
is a small-gap semiconductor/semi-metal, TiS 2 is a semi¬ 
conductor with a bigger gap. 

In conclusion, we have shown that the CDW phase 
in TiSe 2 can be fully predicted by considering only the 
electron-phonon coupling and neglecting any electron 
many-body effect. An incongruence remains between 
the DFT bands and the experiments, especially for the 
high-temperature phase. The simple correction provided 
by the Hubbard-like U term corrects this aspect, but it 
doest not provide a fully coherent picture as it spoils the 
phonon instability. This represents an open issue which 
deserves to be examined in future studies. 
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Appendix A: Folding and Unfolding, definition of 
the spectral weights 

We label with and the direct and reciprocal lattice 
vectors of the unit cell, respectively. Moreover, we define 
& sc and ^ sc as the direct and reciprocal lattice vectors of 
the 2x2x2 supercell, respectively. We find eight vectors 
Gi £ ^sc (defined up to a vector of &) whose differences 
are not vectors of 


Gi £ ^ sc 

Gi -Gj£M 


(Al) 


In the coordinates relative to it is, for example: 


G 0 = (0.0,0.0,0.0) 
G 2 = (0.5,0.0,0.0) 
G 4 = (0.0,0.5,0.0) 
G 6 = (0.5,0.5,0.0) 


Gi = (0.0,0.0,0.5) 
G 3 = (0.5,0.0,0.5) 
G 5 = (0.0,0.5,0.5) 
G 7 = (0.5,0.5,0.5) 


(A2) 


By using these vectors we can unfold a general K £ SBZ 
to eight ki £ BZ: 


ki = K + Gi 


(A3) 
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and write: 

8 

^K,J = Y J V ki (^K,j) (A4) 


In order to plot the unfolded energy spectrum along 
an high-symmetry line of BZ, for each feo of this line 
we considered the superstructure eigenvectors T K(k 0 ),j 
and eigenvalues E K ( ko ^j, where K(ko) is the SBZ point 
where feo G BZ fold into: 


where V ki is the projector on the space of the Bloch 
functions having pseudo-momentum ki with respect to 
that is Vki k,j ) is a function '*/V’ J ( a 0 such that: 

^k i ’ J ( x + r ) = eiki ' r ( x ) (A5) 

The square modulus of V ki k,j ) is, by definition, the 
unfolded weight cof the superlattice band (K, J ) in 
ki : 

<’ J =II^(^)I| 2 (A6) 


K(k 0 ) = ko-G 0 = k 0 (A7) 


Notice that, properly speaking, feo and K(ko) are not 
simple vectors: they represent classes of vectors defined 
up to the sum with an element in 8% and ^ sc , respec¬ 
tively; so even if the representatives feo and K(ko) are 
equal they refer to different set of vectors. Then we plot¬ 
ted the bands Ejc( ko ) i j with intensity: 


r K(k 0 ),J _ , K(k 0 ),J 
J k 0 ~ U ko 


(A8) 


I k ^ k °^ J = 0 corresponding to full transparency (no 
band) and = 1 to full opacity. 
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Supplementary material: 

Electronic and vibrational properties of TiSe 2 in the 
charge-density wave phase from first principles 


GEOMETRICAL PARAMETERS IN THE CHARGE-DENSITY WAVE PHASE 

In this Supplementary Material we report the geometrical parameters for the distorted phase obtained in the studied 
cases. The crystal axes a, c of the hexagonal 2x2x2 supercell are in cartesian coordinates (Angstrom units). The 
atomic positions ai are given in crystal coordinates. For each case, we report both the set of parameters obtained 
first by minimizing the energy along the distortion path (a, c, a 1 ) and then relaxing the structure (a r i x , c r i x , a^ lx ). 


TABLE I. Hexagonal 2 x 2 x 2 supercell parameters in cartesian coordinates for the minimum of the energy along the distortion 
path (a, c) and subsequent relaxed structure (a r i x ,c r i x ). 



a (A) 

C (A) 

«rlx (A) 

Crlx (A) 

LDAexp 

7.080 

12.014 

7.080 

12.014 

gga Bxp 

7.080 

12.014 

7.080 

12.014 

GGAL d P w 

7.080 

12.014 

7.080 

12.014 

GGAxh 

7.072 

13.438 

7.076 

13.456 

GGA™ W 

7.019 

12.331 

7.023 

12.339 


TABLE II. LDAexp atomic positions corresponding to the minimum 
structure (cdix)* The atomic positions are in crystal coordinates. 

along the distortion path (cd) and subsequent relaxed 


a 1 

a 2 

a 3 

a llx 

CLix 

^?ix 

Ti 

-0.0000000 

-0.0085192 

0.0000000 

-0.0000000 

-0.0085705 

0.0000000 

Se 

0.1695741 

0.3362407 

0.1247672 

0.1694345 

0.3362317 

0.1249038 

Se 

0.3333333 

0.1666666 

-0.1247672 

0.3333333 

0.1666666 

-0.1247564 

Ti 

0.0000000 

0.0085192 

0.5000000 

0.0000000 

0.0085705 

0.5000000 

Se 

0.1637590 

0.3304256 

0.6247672 

0.1637681 

0.3305653 

0.6249038 

Se 

0.3333333 

0.1666666 

0.3752328 

0.3333333 

0.1666666 

0.3752436 

Ti 

-0.0000001 

0.4999998 

0.0000000 

-0.0000001 

0.4999998 

0.0000000 

Se 

0.1637592 

0.8333335 

0.1247672 

0.1637683 

0.8332029 

0.1249038 

Se 

0.3333335 

0.6637594 

-0.1247672 

0.3332028 

0.6637684 

-0.1249038 

Ti 

-0.0000001 

0.4999998 

0.5000000 

-0.0000001 

0.4999998 

0.5000000 

Se 

0.1695743 

0.8333335 

0.6247672 

0.1694347 

0.8332029 

0.6249038 

Se 

0.3333335 

0.6695745 

0.3752328 

0.3332028 

0.6694348 

0.3750962 

Ti 

0.5085192 

0.0085192 

0.0000000 

0.5085703 

0.0085704 

-0.0000000 

Se 

0.6666666 

0.3304256 

0.1247672 

0.6667972 

0.3305653 

0.1249038 

Se 

0.8304258 

0.1666666 

-0.1247672 

0.8305654 

0.1667973 

-0.1249038 

Ti 

0.4914808 

-0.0085192 

0.5000000 

0.4914297 

-0.0085704 

0.5000000 

Se 

0.6666666 

0.3362407 

0.6247672 

0.6667972 

0.3362317 

0.6249038 

Se 

0.8362408 

0.1666666 

0.3752328 

0.8362318 

0.1667973 

0.3750962 

Ti 

0.4914806 

0.4999998 

0.0000000 

0.4914296 

0.4999998 

0.0000000 

Se 

0.6666668 

0.8333335 

0.1247672 

0.6666667 

0.8333335 

0.1247564 

Se 

0.8362410 

0.6695745 

-0.1247672 

0.8362320 

0.6694348 

-0.1249038 

Ti 

0.5085191 

0.4999998 

0.5000000 

0.5085702 

0.4999997 

0.5000000 

Se 

0.6666668 

0.8333335 

0.6247672 

0.6666667 

0.8333335 

0.6247564 

Se 

0.8304259 

0.6637594 

0.3752328 

0.8305656 

0.6637684 

0.3750962 
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TABLE III. GGAexp atomic positions corresponding to the minimum along the distortion path ( a z ) and subsequent relaxed 
structure (a£i x ). The atomic positions are in crystal coordinates. 



a 1 

a 2 

a 3 

«rlx 

a ?ix 

C^rlx 

Ti 

-0.0077212 

-0.0077212 

0.0000000 

-0.0077862 

-0.0077862 

0.0000000 

Se 

0.1666666 

0.3363905 

0.1276883 

0.1667914 

0.3363670 

0.1278971 

Se 

0.3363906 

0.1666666 

-0.1276883 

0.3363671 

0.1667914 

-0.1278971 

Ti 

0.0077212 

0.0077212 

0.5000000 

0.0077862 

0.0077862 

0.5000000 

Se 

0.1666666 

0.3302759 

0.6276883 

0.1667914 

0.3304242 

0.6278971 

Se 

0.3302760 

0.1666666 

0.3723117 

0.3304244 

0.1667914 

0.3721029 

Ti 

0.0077211 

0.4999998 

0.0000000 

0.0077861 

0.4999998 

0.0000000 

Se 

0.1666668 

0.8333335 

0.1276883 

0.1666668 

0.8333335 

0.1277357 

Se 

0.3302762 

0.6636096 

-0.1276883 

0.3304245 

0.6636332 

-0.1278971 

Ti 

-0.0077213 

0.4999998 

0.5000000 

-0.0077863 

0.4999998 

0.5000000 

Se 

0.1666668 

0.8333335 

0.6276883 

0.1666668 

0.8333335 

0.6277357 

Se 

0.3363908 

0.6697242 

0.3723117 

0.3363672 

0.6695759 

0.3721029 

Ti 

0.5000000 

0.0077212 

0.0000000 

0.5000000 

0.0077862 

-0.0000000 

Se 

0.6636093 

0.3302759 

0.1276883 

0.6636328 

0.3304242 

0.1278971 

Se 

0.8333333 

0.1666666 

-0.1276883 

0.8333333 

0.1666666 

-0.1277357 

Ti 

0.5000000 

-0.0077212 

0.5000000 

0.5000000 

-0.0077862 

0.5000000 

Se 

0.6697239 

0.3363905 

0.6276883 

0.6695756 

0.3363670 

0.6278971 

Se 

0.8333333 

0.1666666 

0.3723117 

0.8333333 

0.1666666 

0.3722643 

Ti 

0.4999999 

0.4999998 

0.0000000 

0.4999999 

0.4999998 

0.0000000 

Se 

0.6697240 

0.8333335 

0.1276883 

0.6695757 

0.8332087 

0.1278971 

Se 

0.8333335 

0.6697242 

-0.1276883 

0.8332086 

0.6695759 

-0.1278971 

Ti 

0.4999999 

0.4999998 

0.5000000 

0.4999999 

0.4999998 

0.5000000 

Se 

0.6636095 

0.8333335 

0.6276883 

0.6636329 

0.8332087 

0.6278971 

Se 

0.8333335 

0.6636096 

0.3723117 

0.8332087 

0.6636331 

0.3721029 


TABLE IV. GGAex^ atomic positions corresponding to the minimum 
structure (a^ix)* The atomic positions are in crystal coordinates. 

along the distortion path (a 1 ) and subsequent relaxed 


a 1 

a 2 

a 3 

cqix 

Q; r ix 

cqix 

Ti 

-0.0074786 

-0.0074786 

0.0000000 

-0.0075626 

-0.0075626 

0.0000000 

Se 

0.1666666 

0.3362550 

0.1275200 

0.1667657 

0.3362053 

0.1276938 

Se 

0.3362551 

0.1666666 

-0.1275200 

0.3362055 

0.1667657 

-0.1276938 

Ti 

0.0074786 

0.0074786 

0.5000000 

0.0075626 

0.0075626 

0.5000000 

Se 

0.1666666 

0.3304114 

0.6275200 

0.1667657 

0.3305602 

0.6276938 

Se 

0.3304115 

0.1666666 

0.3724800 

0.3305604 

0.1667657 

0.3723062 

Ti 

0.0074785 

0.4999998 

0.0000000 

0.0075625 

0.4999998 

0.0000000 

Se 

0.1666668 

0.8333335 

0.1275200 

0.1666668 

0.8333335 

0.1275652 

Se 

0.3304116 

0.6637451 

-0.1275200 

0.3305605 

0.6637948 

-0.1276937 

Ti 

-0.0074787 

0.4999998 

0.5000000 

-0.0075627 

0.4999998 

0.5000000 

Se 

0.1666668 

0.8333335 

0.6275200 

0.1666668 

0.8333335 

0.6275652 

Se 

0.3362553 

0.6695888 

0.3724800 

0.3362056 

0.6694399 

0.3723063 

Ti 

0.5000000 

0.0074786 

0.0000000 

0.5000000 

0.0075626 

-0.0000000 

Se 

0.6637448 

0.3304114 

0.1275200 

0.6637945 

0.3305602 

0.1276937 

Se 

0.8333333 

0.1666666 

-0.1275200 

0.8333333 

0.1666666 

-0.1275652 

Ti 

0.5000000 

-0.0074786 

0.5000000 

0.5000000 

-0.0075626 

0.5000000 

Se 

0.6695884 

0.3362550 

0.6275200 

0.6694396 

0.3362053 

0.6276937 

Se 

0.8333333 

0.1666666 

0.3724800 

0.8333333 

0.1666666 

0.3724348 

Ti 

0.4999999 

0.4999998 

0.0000000 

0.4999999 

0.4999998 

0.0000000 

Se 

0.6695886 

0.8333335 

0.1275200 

0.6694397 

0.8332344 

0.1276938 

Se 

0.8333335 

0.6695888 

-0.1275200 

0.8332343 

0.6694399 

-0.1276938 

Ti 

0.4999999 

0.4999998 

0.5000000 

0.4999999 

0.4999998 

0.5000000 

Se 

0.6637449 

0.8333335 

0.6275200 

0.6637946 

0.8332344 

0.6276938 

Se 

0.8333335 

0.6637451 

0.3724800 

0.8332343 

0.6637948 

0.3723062 
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TABLE V. GGATh atomic positions corresponding to the minimum along the distortion path (cd) and subsequent relaxed 
structure (a£i x ). The atomic positions are in crystal coordinates. 



a 1 

a 2 

a 3 

«rlx 


a rlx 

Ti 

-0.0105614 

-0.0105614 

0.0000000 

-0.0113010 

-0.0113010 

- 0.0000000 

Se 

0.1666666 

0.3376787 

0.1152080 

0.1672981 

0.3376220 

0.1155628 

Se 

0.3376788 

0.1666666 

-0.1152080 

0.3376221 

0.1672982 

-0.1155628 

Ti 

0.0105614 

0.0105614 

0.5000001 

0.0113010 

0.0113010 

0.5000001 

Se 

0.1666666 

0.3289877 

0.6152078 

0.1672982 

0.3296759 

0.6155626 

Se 

0.3289877 

0.1666666 

0.3847921 

0.3296759 

0.1672982 

0.3844373 

Ti 

0.0105614 

0.4999998 

0.0000000 

0.0113009 

0.4999999 

- 0.0000000 

Se 

0.1666666 

0.8333331 

0.1152080 

0.1666665 

0.8333331 

0.1147726 

Se 

0.3289877 

0.6623209 

-0.1152080 

0.3296759 

0.6623776 

-0.1155628 

Ti 

-0.0105615 

0.4999998 

0.5000001 

-0.0113010 

0.4999997 

0.5000001 

Se 

0.1666666 

0.8333331 

0.6152078 

0.1666665 

0.8333331 

0.6147724 

Se 

0.3376787 

0.6710120 

0.3847921 

0.3376221 

0.6703237 

0.3844373 

Ti 

0.4999999 

0.0105614 

0.0000000 

0.5000000 

0.0113009 

0.0000000 

Se 

0.6623210 

0.3289877 

0.1152080 

0.6623777 

0.3296759 

0.1155628 

Se 

0.8333332 

0.1666666 

-0.1152080 

0.8333333 

0.1666665 

-0.1147726 

Ti 

0.4999999 

-0.0105614 

0.5000001 

0.4999998 

-0.0113009 

0.5000001 

Se 

0.6710121 

0.3376787 

0.6152078 

0.6703238 

0.3376221 

0.6155626 

Se 

0.8333332 

0.1666666 

0.3847921 

0.8333333 

0.1666665 

0.3852275 

Ti 

0.4999999 

0.4999998 

0.0000000 

0.4999998 

0.4999998 

- 0.0000000 

Se 

0.6710120 

0.8333331 

0.1152080 

0.6703239 

0.8327015 

0.1155628 

Se 

0.8333331 

0.6710120 

-0.1152080 

0.8327016 

0.6703238 

-0.1155628 

Ti 

0.4999999 

0.4999998 

0.5000001 

0.4999999 

0.4999999 

0.5000001 

Se 

0.6623210 

0.8333331 

0.6152078 

0.6623777 

0.8327015 

0.6155626 

Se 

0.8333332 

0.6623209 

0.3847921 

0.8327016 

0.6623776 

0.3844373 


TABLE VI. GGAxh W atomic positions corresponding to the minimum along the distortion path (cd) and subsequent relaxed 
structure (a^ix)- The atomic positions are in crystal coordinates. 



a 1 

a 2 

a 3 

^ rlx 



Ti 

-0.0080061 

-0.0080061 

0.0000000 

-0.0081278 

-0.0081278 

0.0000000 

Se 

0.1666667 

0.3365859 

0.1259154 

0.1668046 

0.3365521 

0.1260378 

Se 

0.3365859 

0.1666665 

-0.1259154 

0.3365520 

0.1668043 

-0.1260378 

Ti 

0.0080061 

0.0080061 

0.4999999 

0.0081277 

0.0081277 

0.4999999 

Se 

0.1666668 

0.3300811 

0.6259155 

0.1668046 

0.3302527 

0.6260379 

Se 

0.3300810 

0.1666665 

0.3740845 

0.3302527 

0.1668043 

0.3739622 

Ti 

0.0080059 

0.5000000 

0.0000000 

0.0081277 

0.5000000 

- 0.0000000 

Se 

0.1666665 

0.8333335 

0.1259154 

0.1666666 

0.8333334 

0.1258447 

Se 

0.3300808 

0.6634140 

-0.1259154 

0.3302523 

0.6634478 

-0.1260378 

Ti 

-0.0080063 

0.5000000 

0.4999999 

-0.0081281 

0.4999999 

0.4999999 

Se 

0.1666665 

0.8333335 

0.6259155 

0.1666666 

0.8333335 

0.6258448 

Se 

0.3365856 

0.6699188 

0.3740845 

0.3365518 

0.6697473 

0.3739622 

Ti 

0.4999999 

0.0080061 

0.0000000 

0.4999999 

0.0081280 

0.0000000 

Se 

0.6634142 

0.3300811 

0.1259154 

0.6634480 

0.3302526 

0.1260378 

Se 

0.8333333 

0.1666665 

-0.1259154 

0.8333333 

0.1666665 

-0.1258447 

Ti 

0.4999999 

-0.0080061 

0.4999999 

0.4999998 

-0.0081279 

0.4999999 

Se 

0.6699191 

0.3365859 

0.6259155 

0.6697475 

0.3365521 

0.6260379 

Se 

0.8333333 

0.1666665 

0.3740845 

0.8333333 

0.1666666 

0.3741553 

Ti 

0.5000002 

0.5000000 

0.0000000 

0.5000002 

0.5000000 

0.0000000 

Se 

0.6699194 

0.8333335 

0.1259154 

0.6697477 

0.8331956 

0.1260378 

Se 

0.8333331 

0.6699188 

-0.1259154 

0.8331952 

0.6697472 

-0.1260378 

Ti 

0.5000002 

0.5000000 

0.4999999 

0.5000000 

0.4999998 

0.4999999 

Se 

0.6634145 

0.8333335 

0.6259155 

0.6634484 

0.8331957 

0.6260379 

Se 

0.8333331 

0.6634140 

0.3740845 

0.8331953 

0.6634479 

0.3739622 



